Whole-genome SNP genotyping unveils ancestral and recent introgression in wild and domestic goats.

After the domestication of goats around 10,000 years before the present (BP), humans transported goats far beyond the range of their wild ancestor, the bezoar goat. This brought domestic goats into contact with many wild goat species such as ibex and markhor, enabling introgression between domestic and wild goats. To investigate this, while shedding light on the taxonomic status of wild and domestic goats, we analysed genome-wide SNP data of 613 specimens from 14 taxonomic units, including Capra hircus, C. pyrenaica, C. ibex (from Switzerland, Austria, Germany and Slovenia), C. aegagrus aegagrus, C. a. cretica, C. h. dorcas, C. caucasica caucasica, C. c. severtzovi, C. c. cylindricornis, C. falconeri, C. sibirica sibirica, C. s. alaiana and C. nubiana, as well as Oreamnos americanus (mountain goat) as an outgroup. To trace gene flow between domestic and wild goats, we integrated genotype data of local goat breeds from the Alps as well as from countries such as Spain, Greece, Türkiye, Egypt, Sudan, Iran, Russia (Caucasus and Altai) and Pakistan. Our phylogenetic analyses displayed a clear separation between bezoar-type and ibex-type clades with wild goats from the Greek islands of Crete and Youra clustered within domestic goats, confirming their feral origin. Our analyses also revealed gene flow between the lineages of Caucasian tur and domestic goats that most likely occurred before or during early domestication. Within the clade of domestic goats, analyses inferred gene flow between African and Iberian goats. The detected events of introgression were consistent with previous reports and offered interesting insights into the historical relationships among domestic and wild goats.

Genetic analysis of growth curve in Moghani Sheep using Bayesian and restricted maximum likelihood.

This study was conducted to predict the genetic (co)variance components of growth curve parameters of Moghani sheep breed using the following information: birth weight (N = 7278), 3-mo-old weight (N = 5881), 6-mo-old weight (N = 5013), 9-mo-old weigh (N = 2819], and 12-mo-old weight (N = 2883). The growth parameters (A: maturity weight, B: growth rate, and K: maturity rate) were calculated using Gompertz, Logistic, Brody, and Von Bertalanffy nonlinear models via NLIN procedure of SAS software. The aforementioned models were compared using Akaike information criterion, root mean square error, adjusted co-efficient of determination. Also, both Bayesian (using MTGSAM) and RMEL (using WOMBAT) paradigms were adapted to predict the genetic (co)variance components of growth parameters (A, B, K) due to the best fitted growth models. It was turned out that Von Bertalanffy best fitted to the data in this study. The year of birth and lamb gender had a significant effect on maturity rate (P < 0.01). Also it turned out that within the growth parameter, with increasing (co)variance matrix complexity, the Bayesian paradigm fitted well to the data than the restricted maximum likelihood (REML) one. However, for simple animal model and across all growth parameters, REML outperformed Bayesian. In this way, the h2a predicted (0.15 ± 0.05), (0.11±.05), and (0.04 ± 0.03) for A, B, and K parameters, respectively. Practically, in terms of breeding plan, we could see that genetic improvement of growth parameters in this study is not a tractable strategy to follow up and improvement of the management and environment should be thoroughly considered. In terms of paradigm comparison, REML's bias correction bears up an advantageous approach as far as we are concerned with small sample size. To this end, REML predictions are fairly accurate but the mode of posterior distributions could be overestimated. Finally, the differences between REML and Bayesian estimates were found for all parameter data in this study. We conclude that simulation studies are necessary in order to trade off these parading in the complex random effects scenarios of genetic individual model.

The Iran plateau is known to be the origin of many sheep species nowadays. In Iran, different production systems are operated ranging from intensive to lower-input/extensive ones. However, the majority of these sheep breeds are extensively managed, where lambs are born outside and with little intervention and generally they experience frequent drought and shortage of nutritional value of forages. Meanwhile, the weight of lambs as a whole play a major decision role in rearing or culling them. However, investigations involving the possible genetic improvement of lamb weights over different periods of time have found low genetic variations. This study serves to be comprehensive in addressing this issue in Moghani sheep breed. Fitting many different genetics models over both restricted maximum likelihood and Bayesian paradigms indicated that heritability of weight spanned 0.03 to 0.23. The low genetic variation would lead to recommendations that the improvement of Moghani lamb weights should rather be based upon modification of the environment to create conditions suitable for weight of lambs. This reflects that breeders of Moghani sheep breed have less options to tackle Iran harsh conditions using Moghani sheep genetic potentials.

Whole-Genome Sequencing Data Reveal New Loci Affecting Milk Production in German Black Pied Cattle (DSN).

German Black Pied (DSN) is considered an ancestral population of the Holstein breed. The goal of the current study was to fine-map genomic loci for milk production traits and to provide sequence variants for selection. We studied genome-wide associations for milk-production traits in 2160 DSN cows. Using 11.7 million variants from whole-genome sequencing of 304 representative DSN cattle, we identified 1980 associated variants (-log10(p) ≥ 7.1) in 13 genomic loci on 9 chromosomes. The highest significance was found for the MGST1 region affecting milk fat content (-log10(p) = 11.93, MAF = 0.23, substitution effect of the minor allele (ßMA) = -0.151%). Different from Holstein, DGAT1 was fixed (0.97) for the alanine protein variant for high milk and protein yield. A key gene affecting protein content was CSN1S1 (-log10(p) = 8.47, MAF = 049, ßMA = -0.055%) and the GNG2 region (-log10(p) = 10.48, MAF = 0.34, ßMA = 0.054%). Additionally, we suggest the importance of FGF12 for protein and fat yield, HTR3C for milk yield, TLE4 for milk and protein yield, and TNKS for milk and fat yield. Selection for favored alleles can improve milk yield and composition. With respect to maintaining the dual-purpose type of DSN, unfavored linkage to genes affecting muscularity has to be investigated carefully, before the milk-associated variants can be applied for selection in the small population.

Review: Genetic and protein variants of milk caseins in goats.

The milk casein genes in goats, are highly polymorphic genes with numerous synonymous and non-synonymous mutations. So far, 20 protein variants have been reported in goats for alpha-S1-casein, eight for beta-casein, 14 for alpha-S2-casein, and 24 for kappa-casein. This review provides a comprehensive overview on identified milk casein protein variants in goat and non-coding DNA sequence variants with some affecting the expression of the casein genes. The high frequency of some casein protein variants in different goat breeds and geographical regions might reflect specific breeding goals with respect to milk processing characteristics, properties for human nutrition and health, or adaptation to the environment. Because protein names, alongside the discovery of protein variants, go through a historical process, we linked old protein names with new ones that reveal more genetic variability. The haplotypes across the cluster of the four genetically linked casein genes are recommended as a valuable genetic tool for discrimination between breeds, managing genetic diversity within and between goat populations, and breeding strategies. The enormous variation in the casein proteins and genes is crucial for producing milk and dairy products with different properties for human health and nutrition, and for genetic improvement depending on local breeding goals.

Validation of somatic cell score-associated SNPs from Holstein cattle in Sudanese Butana and Butana × Holstein crossbred cattle.

The Bos indicus zebu cattle Butana is the most commonly used indigenous dairy cattle breed in Sudan. In the last years, high-yielding Holstein dairy cattle were introgressed into Butana cattle to improve their milk yield and simultaneously keep their good adaption to extreme environmental conditions. With the focus on the improvement of milk production, other problems arose such as an increased susceptibility to mastitis. Thus, genetic selection for mastitis resistance should be considered to maintain healthy and productive cows. In this study, we tested 10 single nucleotide polymorphisms (SNPs) which had been associated with somatic cell score (SCS) in Holstein cattle for association with SCS in 37 purebred Butana and 203 Butana × Holstein crossbred cattle from Sudan. Animals were genotyped by competitive allele-specific PCR assays and association analysis was performed using a linear mixed model. All 10 SNPs were segregating in the crossbred Butana × Holstein populations, but only 8 SNPs in Sudanese purebred Butana cattle. The SNP on chromosome 13 was suggestively associated with SCS in the Butana × Holstein crossbred population (rs109441194, 13:79,365,467, PBF = 0.054) and the SNP on chromosome 19 was significantly associated with SCS in both populations (rs41257403, 19:50,027,458, Butana: PBF = 0.003, Butana × Holstein: PBF = 6.2 × 10-16). The minor allele of both SNPs showed an increase in SCS. Therefore, selection against the disadvantageous minor allele could be used for genetic improvement of mastitis resistance in the studied populations. However, investigations in a bigger population and across the whole genome are needed to identify additional genomic loci.

Capture Sequencing to Explore and Map Rare Casein Variants in Goats.

Genetic variations in the four casein genes CSN1S1, CSN2, CSN1S2, and CSN3 have obtained substantial attention since they affect the milk protein yield, milk composition, cheese processing properties, and digestibility as well as tolerance in human nutrition. Furthermore, milk protein variants are used for breed characterization, biodiversity, and phylogenetic studies. The current study aimed at the identification of casein protein variants in five domestic goat breeds from Sudan (Nubian, Desert, Nilotic, Taggar, and Saanen) and three wild goat species [Capra aegagrus aegagrus (Bezoar ibex), Capra nubiana (Nubian ibex), and Capra ibex (Alpine ibex)]. High-density capture sequencing of 33 goats identified in total 22 non-synonymous and 13 synonymous single nucleotide polymorphisms (SNPs), of which nine non-synonymous and seven synonymous SNPs are new. In the CSN1S1 gene, the new non-synonymous SNP ss7213522403 segregated in Alpine ibex. In the CSN2 gene, the new non-synonymous SNPs ss7213522526, ss7213522558, and ss7213522487 were found exclusively in Nubian and Alpine ibex. In the CSN1S2 gene, the new non-synonymous SNPs ss7213522477, ss7213522549, and ss7213522575 were found in Nubian ibex only. In the CSN3 gene, the non-synonymous SNPs ss7213522604 and ss7213522610 were found in Alpine ibex. The identified DNA sequence variants led to the detection of nine new casein protein variants. New variants were detected for alpha S1 casein in Saanen goats (CSN1S1 ∗C1), Bezoar ibex (CSN1S1 ∗J), and Alpine ibex (CSN1S1 ∗K), for beta and kappa caseins in Alpine ibex (CSN2 ∗F and CSN3 ∗X), and for alpha S2 casein in all domesticated and wild goats (CSN1S2 ∗H), in Nubian and Desert goats (CSN1S2 ∗I), or in Nubian ibex only (CSN1S2 ∗J and CSN1S2 ∗K). The results show that most novel SNPs and protein variants occur in the critically endangered Nubian ibex. This highlights the importance of the preservation of this endangered breed. Furthermore, we suggest validating and further characterizing the new casein protein variants.

Whey protein polymorphisms in Sudanese goat breeds.

The aim of the present study was to assess genetic variation that is characteristic for Sudanese goat breeds in the milk whey protein genes (LALBA and BLG). Four Sudanese goat breeds were screened for variability in LALBA and BLG genes at the DNA level by comparative sequencing of five animals per breed. Sixteen SNPs were identified in LALBA: seven in the upstream region, six synonymous, and three in the 3´-UTR. Three novel synonymous SNPs in exon 2 (ss5197800003, ss5197800012, and ss5197800004) were found in Nubian, Desert, and Nilotic, but not in Taggar goats. One SNP in the promoter of LALBA (rs642745519) modifies a predicted transcription factor binding site for Tcfe2a. The SNPs in the 3'-UTR (rs657915405, rs641559728, and rs664225585) affect predicted miRNA target sites. With respect to haplotypes in the exonic region, haplotype LALBA-A is most frequent in Nubian, Desert, and Nilotic goats, while haplotype LALBA-D is prevalent in Taggar goats. In BLG, 30 SNPs were detected: eight in the upstream gene region, two synonymous, 17 intronic, and three in the 3'-UTR. Among the 30 identified SNPs, 15 were novel. Four of these novel SNPs were located in the upstream gene region, one was synonymous, and ten were intronic. The novel synonymous SNP (ss5197800017), located in exon 2, was only found in Nubian and Nilotic goats. The SNPs ss5197800010 and rs635615192 in the promoter are located in predicted binding sites of transcription factors (M6097, Elk3, Elf5, and GABPA). Among seven haplotypes detected in the coding region, haplotype BLG-A is most frequent in Nubian and Nilotic goats while haplotype BLG-B is most frequent in Desert and Taggar goats. The high variability in regulatory gene regions among Sudanese goats could potentially affect the quality and yield of whey proteins in goat milk and provide a wide resource for genetic improvement of milk production and milk technology characteristics.

Whole genome population genetics analysis of Sudanese goats identifies regions harboring genes associated with major traits.

BACKGROUND: Sudan is endowed with a variety of indigenous goat breeds which are used for meat and milk production and which are well adapted to the local environment. The aim of the present study was to determine the genetic diversity and relationship within and between the four main Sudanese breeds of Nubian, Desert, Taggar and Nilotic goats. Using the 50 K SNP chip, 24 animals of each breed were genotyped. RESULTS: More than 96% of high quality SNPs were polymorphic with an average minor allele frequency of 0.3. In all breeds, no significant difference between observed (0.4) and expected (0.4) heterozygosity was found and the inbreeding coefficients (FIS) did not differ from zero. Fst coefficients for the genetic distance between breeds also did not significantly deviate from zero. In addition, the analysis of molecular variance revealed that 93% of the total variance in the examined population can be explained by differences among individuals, while only 7% result from differences between the breeds. These findings provide evidence for high genetic diversity and little inbreeding within breeds on one hand, and low diversity between breeds on the other hand. Further examinations using Nei's genetic distance and STRUCTURE analysis clustered Taggar goats distinct from the other breeds. In a principal component (PC) analysis, PC1 could separate Taggar, Nilotic and a mix of Nubian and Desert goats into three groups. The SNPs that contributed strongly to PC1 showed high Fst values in Taggar goat versus the other goat breeds. PCA allowed us to identify target genomic regions which contain genes known to influence growth, development, bone formation and the immune system. CONCLUSIONS: The information on the genetic variability and diversity in this study confirmed that Taggar goat is genetically different from the other goat breeds in Sudan. The SNPs identified by the first principal components show high Fst values in Taggar goat and allowed to identify candidate genes which can be used in the development of breed selection programs to improve local breeds and find genetic factors contributing to the adaptation to harsh environments.

The F279Y polymorphism of the GHR gene and its relation to milk production and somatic cell score in German Holstein dairy cattle.

The bovine growth hormone receptor (GHR) gene has been identified as a strong positional and functional candidate gene influencing milk production. A non-synonymous single nucleotide polymorphism (SNP) in exon 8 leads to a phenylalanine to tyrosine amino acid substitution (F279Y) in the receptor. The aim of the study was to estimate the effects of the F279Y mutation on milk yield, fat, protein, casein, and lactose yield and content, as well as somatic cell score (SCS), in a German Holstein dairy cattle population. The analysis of 1,370 dairy cows confirmed a strong association of the F279Y polymorphism with milk yield, as well as with fat, protein, and casein contents. Furthermore, increasing effects on lactose yield and content for the 279Y allele were found. Even though the tyrosine variant occurred as the minor allele (16.5%), its substitution effects were 320 kg (305 d), 0.02 kg per day, and 0.07 kg per day for milk, casein, and lactose yields, respectively. The same allele had negative effects on fat, protein, and casein contents. Finally, the high-milk-yield tyrosine allele was also associated with lower SCS (p < 0.05). The data support the high potential of the F279Y polymorphism as a marker for the improvement of milk traits in selection programs.